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It is shown how a quantum cellular automaton can describe very precisely the Dirac evolution, 
without requiring Lorentz covariance. The automaton is derived with the only assumptions of 
minimal dimension and parity and time-reversal invariance. The automaton extends the Dirac field 
theory to the Planck and ultrarelativistic scales. The Dirac equation is recovered in the usual 
particle physics scale of inertial mass and momenta. In this first paper the simplest case of one 
space dimension is analyzed. We provide a technique to derive an analytical approximation of 
the evolution of the automaton in terms of a momentum-dependent Schrodinger equation. Such 
approximation works very well in all regimes, including ultrarelativistic and Planckian, for the 
typical smooth quantum states of field theory with limited bandwidth in momentum. Finally we 
discuss some thought experiments for falsifying the existence of the automaton at the Planck scale. 
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I. INTRODUCTION 

Quantum Field Theory (QFT) is the most detailed de- 
scription of the dynamics of physical systems available. 
Its predictions are extremely accurate, and have been 
tested in a countless number of experiments. Despite its 
impressive predictive power, QFT is affected by a num- 
ber of conceptually relevant issues, including all problems 
arising from the continuum, localizability, causality vio- 
lation, and quantization (for a review, see Refs. [IH5]). 
The main framework is made of recipes that are validated 
a posteriori, e.g. the quantization rules, both the canon- 
ical and the Feynman's path integral, the former being 
affected by the ordering issue, the latter being based on 
a mathematically undefined notion. In practice a quan- 
tized theory is first formulated, and then checked a pos- 
teriori if it is renormalizable, otherwise it is discarded. 
This process has been one of the main theory-sieve in 
the formulation of the so-called Standard Model [§], the 
worldwide collaborative effort that succeeded in explain- 
ing a wide variety of experimental results, to the extent 
of being regarded as a "theory of almost everything". 
However, the price to pay is the limited validity of the 
theory, which does not accomodate gravity. 

The issue of localizability [TJ [7J in QFT-the essential 
ingredient of the particle concept-is widely considered 
one of the most severe, since it prevents the formula- 
tion of a rigorous quantum theory of measurement for 
fields. It is also the main hallmark of the unsolvable 
tension between relativity and quantum theory, as for 
e.g. the boost-induced delocalization of Newton- Wigner 
states [8|. This is also another side of the causality vi- 
olation due to the Hamiltonian description, which leads 
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to the wave- function superluminal tails [5]. 

There is no easy way for curing all the issues plaguing 
QFT, if not through a reconsideration of its very founda- 
tions. It is common opinion that the continuum descrip- 
tion of space-time is the source of the problem, and that 
at some very small scale-the Planck scale-the description 
should be taken discrete. This is the mainstream philos- 
ophy of the approaches to quantum space-time jTOl [H] , 
of loop quantum-gravity [j^MIl] , and, in some way, of the 
causal set approach [15 . The discreteness of the Planck- 
scale has also the consequence that Lorentz covariance 
and all symmetries of QFT are no longer valid, and are 
recovered at the Fermi scale [TorfTS] . Recently, potential 
experimental tests of Planck-scale phenomenology have 
been proposed by several authors [I9H22] . 

A description of the quantum field as a quantum cel- 
lular automaton (QCA) represents a logically coherent 
way of tackling the Planck scale discreteness with min- 
imal assumptions. Compared to other approaches, the 
QCA is motivated by a list of reasons. First: it relaxes 
the tension between quantum theory and relativity, since 
it does not need relativity, which in turn is emergent from 
the former. Second: it is quantum ab-initio, without the 
need of quantization rules. Third: it is free from all the 
problems arising from the continuum, and doesn't suffer 
violations of causality. 

A first classical automaton description for the Dirac 
field evolution is hidden in the nonstandard analysis for- 
mulation of the path-integral of Nakamura [53] , where the 
"infinitesimal" plays the role of the Planck length. Then, 
the word "automaton" first appeared in relation to rela- 
tivistic field-theory in the pioneering work of Bialynicki- 
Birula [24 , where the automaton describes a discretiza- 
tion of the Weyl differential equation. The automaton 
is synthesized by a unitary matrix, but it is still a clas- 
sical one, representing an updating rule on a lattice of 
classical variables. The possibility of using automata for 
describing the evolution of classical relativistic fields also 
emerged in the context of lattice-gas simulations, espe- 
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cially in the seminal work of Meyer [35] , where a notion 
of "field automaton" first appeared, and in the papers 
of Yepez [33]. In all cases, however, the automata are 
classical ones, and are devised as simulation tools for a 
discretized path-integral approach. 

Recently the quantum cellular automaton has been 
proposed as a minimal-assumption extension to the 
Planck and ultrarelativistic scales of QFT [27H30] . A 
Dirac automaton in one-dimension has been derived 
without imposing Lorentz covariance, starting from prin- 
ciples of information-theoretical nature (for a nontechni- 
cal review see Ref.[3T]). As we will see in this paper, the 
automaton is described by a simple 2x2 matrix: how- 
ever, despite its simplicity, the QCA describes a very rich 
physical phenomenology, leading to unexpected interest- 
ing predictions, e.g. it anticipates an upper bound to the 
inertial mass of the Dirac field [2"9"] , as a simple conse- 
quence of unitarity of quantum evolution, without invok- 
ing black-hole arguments from general-relativity. More- 
over, the automaton framework even allows to redefine 
fundamental physical notions and physical constants in 
an informational way |28L I30j . 

The quantum cellular automaton represents an exten- 
sion of QFT, in the sense that it describes also local- 
ized states and measurements that are not manageable 
by QFT, giving a unified description of the field dynam- 
ics at all scales, ranging from the Planck, ultrarelativistic 
to the usual particle-physics scale. 

In this work we provide a thorough study of the 
physics emerging from the Dirac automaton in one space- 
dimension of Ref . |29j , showing how the Dirac dynamics 
is perfectly recovered at the Fermi scale, though rela- 
tivistic covariance and other symmetries are violated at 
the Planck and ultrarelativistic ones (a Dirac quantum 
automaton has also been derived recently for space di- 
mension d — 2,3, and will be the object of the next 
publication [32] of this series). 

The QCA generalizes the notion of cellular automa- 
ton of von Neumann [33] to the quantum case. The 
earliest appearance of the notion of QCA in the liter- 
ature dates back to 1982 from R. Feynman [34 . The 
first QCA as we know it nowadays has been introduced 
in Ref. [35]. In the following literature the QCA has 
been mostly a computer-science object of investigation, 
especially in the field of quantum information, with in- 
teresting general mathematical results [361438] . However, 
the QCA has never been used as a theoretical framework 
for physics, as e.g. field theory, as we will do in this 
paper. 

The classical cellular automaton consists of a regular 
lattice, whose sites can be in a finite number of states 
(e.g. on and off) along with a rule that updates the 
state from "time" t to t + 1 . This rule must be local, 
namely the state at site x at time t + 1 depends only 
on the states of a finite set of neighbouring sites at time 
t. In the quantum version, the sites of the lattice are 
quantum systems interacting unitarily with a finite set of 
neighbouring systems. In the present context, the lattice 



is a Planck scale version of the notion of space. The 
metric corresponds tocounting systems along the lattice, 
and the conversion from counting to the usual metric is 
given by a unit length which is equal to the Planck length 
£p. In a similar way the counting of the updating steps 
corresponds to time, with the conversion factor given by 
the Planck time rp. Specifically as a framework for QFT 
the quantum systems are field evaluations on the lattice, 
and a priori it can be taken Bosonic or Fermionic. 

After reviewing the one-dimensional Dirac automaton, 
in Sect. [IT] we show how the automaton recovers pre- 
cisely the Dirac dynamics in the limit of small masses 



and momenta. Then in Sect III we present an analyti- 
cal approximation method for evaluating the automaton 
evolution for single-particle states smooth in momentum 
and with limited bandwidth-the typical quantum states 
of field theory. In this way we derive a momentum- 
dependent Schrodinger equation, which works very well 
in all regimes, including ultrarelativistic and Planckian. 
We compare computer simulations with the analytic ap- 
proximation, and provide the leading order corrections 
to the Dirac equation. After discussing thought exper- 
iments for falsifying the existence of the automaton at 



the Planck scale in Sect IV we conclude the paper with 
future perspectives. 



II. THE ONE-DIMENSIONAL DIRAC 
AUTOMATON 

In this section we review the one-dimensional Dirac 
automaton of Ref. [55] , which is the focus of the present 
paper. Then we show how the Dirac dynamics is recov- 
ered for low momenta and masses. 



A. Derivation of the Dirac Automaton 

In QFT the unitary evolution of a field ip(x) with A 
internal degrees of freedom is given by 

^{t + l) = uH(t)U, il>{x) = {V>a(x)}a=i,...,a. (1) 

The QCA corresponding to the one-step update of the 
field can be described by a unitary matrix U acting lin- 
early on the field operator i/> 



V>(t + i) = u^(t), V(<) 



i/>(x,t) 
ip(x + l,t) 



(2) 



where we use the symbols x £ Z and t E Z for the adi- 
mensional numbering of the lattice sites and time steps 
respectively. A pictorial representation of the automaton 
U is shown in Fig jl] where the light cones depict the local 
structure of the evolution, and at every site x € Z of the 
lattice there are A systems corresponding to the internal 
components of the field. 
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When the lattice dimension is one, the most general 
translational invariant QCA with next-neighbouring in- 
teraction is given by 



U = R5 + L5 f + M, 

where S denotes the shift operator 

Si>(x) := tp(x + l), 



(3) 



(4) 



and R, L and M are A x A matrixes which do not depend 
on the lattice site, due to translation invariance. The 
unitarity of U implies 



MR 



RR f 

hLM f 



LL T + MM* = I 



0. 



LR 







(5) 



We now derive the Dirac automaton as the minimal- 
dimension QCA, satisfying the following requirements 

i. Unitarity of the evolution; 

ii. Homogeneity of the interaction topology; 

iii. Invariance under time-reversal 1 1— > —t; 

iv. Invariance under parity x i— > — x; 

v. Minimal dimension for a non-identical evolution. 

The first two assumptions are already contained in the 
definition itself of QCA. Assumptions [m] and [Iv] are the 
symmetries of the pure topology of the causal network 
(in Fig. [2]), corresponding to considering the network as 
undressed. The next-neighbouring interaction is not an 
assumption by itself, since it is always possible to reduce 
to such a case by grouping a periodic pattern of the net- 
work into a single node of the automaton. Therefore, the 
only assumption is the minimality one |vj In equations, 
assumptions [m] and [ry] correspond to 



TUT 1 = U f , 
PUP f = R5 f + L5 + M, 



(6) 



where T is the anti-unitary operator associated to the 
time reversal transformation, and P is the unitary op- 
erator of the lattice reflection. Notice that PUP^ is the 
same as U with S exchanged with . For A = 1 the only 
translational invariant QCA satisfying parity invariance 
is the identical one U = I. Next, we have the case A = 2. 
Eq. ([7]) shows that R and L are unitarily equivalent, 
whence, from Eq (Jsj) it follows that they are both rank 
one. Thus we can choose the basis where 



R 



a 1 




ip(x) 



(8) 



naming the two components of the field ipR and ipi, right 
and left modes. We now require P and T to be repre- 
sented in such basis as 



1 

1 



T = C 



1 

1 



(9) 




t + 1 



t 1 r 

tl>(x - 2) tj>{x - 1) i>{x) il>(x + 1) i>(x + 2) 



FIG. 1: Illustration of a one-dimensional quantum cellular 
automaton unitary step. Each site of the lattice x corresponds 
to a quantum field evaluation t/}(x). The field operator at site 
x interacts with the field tp(x ± 1) at neighbouring sites. In 
the case of the Dirac automaton the field operator has two 
components (see text). 




FIG. 2: Schematic of the three time steps causal network cor- 
responding to a one-dimensional quantum cellular automaton 
with next-neighbouring interaction. The topology of the net- 
work is left invariant by the mappings t B -< and x i->- — x 
and the dynamics of the automaton is assumed to be time 
reversal (T) and parity (P) invariant (see Eqs. |6| and Q). 



where C is the anti-unitary operator denoting complex 
conjugation in the representation (|8|. With this choice of 
basis the parity invariance Q, which implies PRP^ = L 
and PMP t = M, 



gives 



R 



ai a 2 






a 2 ai 



M 



a 3 a 4 
04 a 3 



(10) 



with a.i 6 C for i = 1, 



, 4. From the time time-reversal 



invariance that is TRT f = L f and TMT f = M f , it 



follows a 2 = a 3 = 0. Finally, using the unitarity of U ^ 
we get |ai| 2 + |a4 1 2 = 1 and ^(0104) = which, up to a 
global phase, give the unique automaton 



U = 



iS 



m z = 1 



(11) 



The constants n and m in the last equation can be cho- 
sen positive (the relative sign corresponds to take the 
left-mode redefined with a mi nus sign) . The unitarity 
constraint n 2 + m 2 = 1 in Eq. ( 11 ) forces the parameter 
m to be m £ [0, 1]. As we will see in the following, m 
plays the role of an adimensional inertial mass, whereas 
n is the inverse of a refraction index of vacuum. Thus, 
the inertial mass of the field is bounded from above [29] , 
as mentioned in the introduction. In the digital-analog 
conversion from the automaton to the usual dimensional 
Dirac equation, we take the Planck length lp and the 
Planck time rp as conversion factors for space and time 
respectively, whereas the mass conversion factor is taken 
equal to the Planck mass uip. The maximal speed of 
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FIG. 3: (Colors online) The automaton (m = 0.92) evolu- 
tion for a localized state (top) and a smooth state (bottom). 
The localized state cannot be described by QFT, whereas the 
smooth state is a typical QFT state (for more details see text). 
The localized state is given by 1 30) <g> ;4=(1, 1). The smooth 

state is given by ^2 X g(x)e tkx \x) Cg> \+}h with g(x) Gaussian 
with mean value xq — 30 and width a — 3. |+)fc is a one par- 



ticle eigenstate (21 1 of the Dirac automaton with momentum 
k — 0.3-7T (see text). 



local-state propagation in the dimensional case is then 
c = £ p /t p , corresponding to the speed of light. £ P , r P , 
and mp are tree fundamental universal constants defin- 
ing the measure for dimensions [L] [T] [M], and from 
them one can derive any other universal constant, e.g. 
h = rripipc [30] . The automaton for m — corresponds 
to the Weyl equation, and will be refereed as Weyl au- 
tomaton. 

In Fig. [3] we give two computer evaluations of the au- 
tomaton evolution for a localized state not describable by 
QFT, and for a smooth state of the kind used in QFT (for 
more details see the following). The same Dirac automa- 
ton of ( 11 ) is also derived in Ref. [32], by recovering the 

3 as the unique minimal QCA, 



Weyl automaton for d 
whereas the Dirac automaton is then the only possible 



local coupling of two Weyl automata (the only possible 
automata that can be coupled are a pair of reciprocally 
inverse automata). In this way one obtains a two-spinor 
automaton, which for d — 1 splits into two identical au- 
tomata of the form given in Eq (11 ). Originally the au- 



tomaton of (111 has been derived heuristically as the one 



describing the free flow of quantum information |29j . 



B. Recovering the Dirac dynamics 



Here we show how the Dirac QCA ( 11 ) recovers, for low 



momenta and masses, the dynamics of the Dirac equation 
which is 



ihd t ip(x, t) 



-ihcd x 
mc 2 



Zc8 x J**'*)- ( 12 ) 



The digital version of the equation corresponds to us- 
ing adimensional x, t, and m (corresponding to Planck 
units), with h = c = 1. For convenience, as usual in the 
literature (EMEU), we study the dynamics of the au- 
tomaton in the momentum representation. For the field 
operator we have 



1 



/2tt 



Y,e- ikx ^{x), fce[-7r,7r], (13) 



x£2 



where with little abuse of notation we utilize the variable 
name k to denote the Fourier transform of any function of 
x. Notice that the automaton model is naturally band- 
limited k g [— 7r,7r] (k G [— ir/£p, n/£p] in usual units) 
and periodic in momenta due to the discreteness of the 
lattice. Correspondingly, for the unitary matrix U in Eq. 



(Ill we have 



U 



dkU(k), U(jfe) 



ne lk 

—im ne 



nn 

— ik 



(14) 



By taking a real power of the unitary matrix, we define an 
abstract Hamiltonian that describes the automaton evo- 
lution for continuous times, interpolating between time- 
steps, namely U* = exp(— iHt). Upon diagonalizing the 



matrix U(fc) in Eq. (14), one obtains 

H= I dkU(k), 

Tj/ju\ _ w /— nsin(fc) m 
~sin(w) V 171 rasm(fc) 



(15) 
(16) 



where ui(k, m) is the dispersion relation of the automaton 

Lj(k,m) = arccos(V 1 — m 2 cos(fc)). (17) 

Notice that the Hamiltonian H is unphysical, as is un- 
physical the continuous process between two following 
time-steps. Indeed, one could easily see that the Hamil- 
tonian (15 1 involves interactions between all systems, in- 



cluding those very far apart (the "local" Hamiltonian 
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corresponding to the finite-difference form of the usual 
Dirac field Hamiltonian is given by the imaginary part of 

We now expand H(fc) near k = and m = 0, in order 
to recover the Dirac equation for the quantum field, along 
with the first corrections. One has 



H(fc) 



OO 

y - 



i,3=0 



SH(fc) 

dk i dmJ> 



k 1 



(18) 



k,m—0 



and collecting all terms with the same overall order in m 
and k, we find 



H(fc) = H D (fc) 



mk \(k 2 + m 2 )\ 
3 \^(k 2 +m 2 ) -mk J (19) 

0(m p k q ) where p + q = 5 



where Hrj(fc) is the adimensional Dirac Hamiltonian of 
Eq. ( 12 ) in the momentum representation 



H D (fc) = 



— k m 
m k 



(20) 



The Hamiltonian ( 20 1 identifies the parameters k and 
?7i of the automaton with the momentum and the mass 
of the Dirac field, respectively. With the automaton 
operating at the Planck scale the maximal mass must 
coincide with the Planck mass mp. Remembering the 
value nip in usual units mp = 2.17651(13) x 10 _8 kg 
= 1.2209 x 10 28 eV/c 2 , a particle having rest mass 100 
GeV/c 2 as e.g. the Higg's boson corresponds to an adi- 
mensional mass m « 10~ 17 (for the electron and the 
proton one would have m e « 10~ 23 and m p « 10 -19 , 
respectively). Similarly a momentum of 10 20 eV/c, the 
highest ever observed for Ultra-High Energy Cosmic Rays 
(UHECRs, see for example [44]), corresponds in Planck- 
ian units to kcp. ~ 10~ 8 (a high energy momentum of 
ITev in reachable in a LHC experiments corresponds to 
khHC ~ 10~ 16 ). Therefore, the relativistic approxima- 
tion works very well in particle physics experiments. 



C. Dispersion relation 

The eigenvalues and the eigenvectors of the unitary 
matrix U(fc) in Eq. (14) are given by 



Ufe(s) 



1 

•J2 



v 7 ! — sv 
s-\/l + sv 



s = ±, (21) 



in terms of the automaton dispersion relation (17) and 
the group velocity v — dkOJ. 

In Fig. [4] the automaton dispersion relation ui is com- 
pared with the Dirac one 



UJ D 



:= woffc, m) = \Jk 2 + m 2 . 



(22) 



The dispersion relations show an overlap vs k around 
k = for small m corresponding to the typical particle- 
physics regime. Analytically, the leading term correction 



to the Dirac dispersion is given by 

m 2 k 2 — m 2 



u D 1 



6 k 2 



m A 



(23) 



Notice that the smallest non zero value of the mass 
m — 0.3 reported in Fig. [4] is huge compared to any 
particle mass, and still is possible to find a sector in the 
momentum space around k = where the automaton dis- 
persion relation well overlaps the Dirac one. This means 
that for m sufficiently small it is possible to fix a cutoff 
in the momenta \k\ < k < t: such that in that regime 
the Dirac automaton and the usual Dirac evolution are 
very similar. The zone-border value k = tt describes the 
Planck-scale behaviour, whereas smaller k correspond to 
look at the automaton at larger scales. We emphasize 
the conceptual difference between the large scale limit 
analyzed here \k\ < k < tt with £ — lp, r = tp, and 
the continuum limit £ — > 0, t — > considered in Refs. 

MM- 



D. Particle states 

The QCA automaton U generally operates on the vec- 
tor field tp which describes an arbitrary number of par- 
ticles. The vacuum state for the automaton is defined as 
the state such that 



4> s (k)\n) = v.s = ±, vfc g [-7T, tt] 



(24) 



Up to now, we have not specified the nature 
Fermionic/Bosonic of the field, and indeed the same au- 
tomaton could be used to describe both cases (along with 
anyons and parastatistics): the relation spin-statistics 
will be considered for the the automaton with d = 3 
space-dimensions [32]. Here we will focus only on the 
Fermionic case of anticommuting field. A ./V-particle 
state can be obtained by acting with the field operator 
on the vacuum as follows 



(25) 



Specifically, for N — 1 particle eigenstates of U, we write 
^ s (k)\Q) = | S ) fe |fc), (26) 
whereas for N — 2 we have 



^(WWIfi) = \si)Mk 2 \h,k 2 ) (27) 

where |fci,fc2) = —\k2,kx), and so forth for N > 2. The 
corresponding eigenvalues of the (logarithm of ) U are 



JY 



j(N, k, s) = ^2 s i m )- 



(28) 
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FIG. 4: (Colors online) Comparison between the dispersion 
relations u(k, m) of the Dirac automaton and of the Dirac 
equation, in Eqs \17\ and respectively. In the top fig- 

ure the dispersion relation is plotted versus the adimensional 
mass m G [0, 1] and momentum k £ [— 7T, n] (m = 1 corre- 
sponds to the Planck mass). The green surface represent the 
automaton, whereas the blue the Dirac one. In the bottom 
figures u)(k,m) is plotted versus k for four values of m (the 
red line corresponds to the automaton, whereas the black one 
is the Dirac's). We can see that the two dispersion relations 
coincide for small masses and momenta, and the larger the 
mass the smaller the overlap region around k — 0. 



III. THE QUANTUM-FIELD LIMIT OF THE 
AUTOMATON 



The ultimate aim of this approach to QFT is clearly the 
connection to phenomenology. In the last section we will 
provide a bound from an optimization over all possible 
experimental setups of the probability of discriminating 
between the Dirac automaton and usual Dirac evolution. 
This is a powerful result, both from the theoretical and 
the phenomenological point of view, since it allows to 
verify the theoretical discriminability between the two 
theories, giving insights on the most suitable experimen- 
tal scenario where searching for violations of the usual 
Dirac evolution. However, the huge class of experimen- 
tal setups considered in the derivation of the bound in- 
cludes also configurations non realizable with the present 
technology. 

In this section we explore the behaviour of the Dirac 



QCA for the one-particle states of quantum field theory. 
We consider initial states whose momentum distribution 
is smoothly peaked around some k , namely 



Wo)>=/ -= 5 (fc,o)| s ) fc |fc), 



'2tt' 



B=±, 



(29) 



where g(k, 0) € C^°[—tt,tt] is a smooth function satisfy- 
ing the bound 



- / dk\g(k,0)\ 2 > 1-e, <r,e > 0. 

if Jk -a 



(30) 



where the two-component vector Is)/, is defined in Eq. 



(211 



At time t and in the position representation, the state 
in Eq. ((29]) can be written as 



\tp(x,t)) := e l{koX ~ su,ot) \(t){x,t)) 
dk 



i(K x— sQ,(k ,m)t) 



2tt 



g(k,o)\ 



S)k 



where we posed 
K = k — fcoj 

Q(k,m) — uj(k,m) — ljq, ujq = uj(ko,m). 



(31) 



(32) 



A. The fc-dependent Schrodinger equation 

It is convenient to take x, t to be real-valued continuous 



variable by extending the Fourier transform in Eq. (31) 



to real x,t. We derive the integral in Eq. (31) with 
respect to t, and expand fl vs k around ko up to the 
second order. Then, taking the resulting derivatives with 
respect to x out of the integral (using the dominated 
derivative theorem), we obtain the following fc-dependent 
Schrodinger equation with drift 



id t \(f>(x,t)) 



d_ 

dx 



with the drift constant v and the diffusion constant D 
depending on k and m as follows 



1 



D := 



1 + m 2 cot 2 (fc )' 
VT — m 2 m 2 cos (fco) 



(34) 
(35) 



(sin 2 (fc ) +m 2 cos 2 (fco)) 2 
and with the identification of the initial condition 

|^(x,0)) = |^(x,0)). (36) 



The drift and diffusion coefficients arc obtained as deriva- 
tives of the dispersion relation as v = ui^J and D = cu^j , 
where 



(„) _ d n uj{k,m) 

W fco - Q k n 



(37) 
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For \il>(x,Q)) satisfying Eq. (30), Eq. (33) provides the 



approximation of the state of the particle \tp(x,t)) 



O i(k x—suot) 



\(/)(x,t)), corresponding to 



dk 

t 

v/2tt 



-is{ujQ-\-v k— ^Dk 2 )t 



g(k,0)\s) k \k). (38) 



B. Accuracy of the approximation 

The accuracy of this approximation can be quantified 
in terms of the parameters a and e of the initial state 
by evaluating (see Appendix [A]) the overlap between the 
states pil and p8| 



\m)\m)\>i 



<ya 3 t 



0(a 5 )t. (39) 



where 



7 



2vr 



k +a 



dk\g{kM 



kn — c 



Notice that the choice of smooth function g(k,t) cor- 
responds to a two-component vector \ip{x,t)) that is 
Schwartz on the real line, as for the field of QFT. We are 
thus evaluating the asymptotic field-limit of the QCA. 
The approximation (38) works well in all regimes, 



namely for any mass m and momentum k. We can test 
the accuracy of the approximation by comparing \ip{x, t)) 
with the automaton simulation for given initial state. In 
Fig. [5] we show an example in the Planckian ultrarela- 
tivistic regime. Here the Schwartz-class field state (29) 



is a superposition of Hcrmite functions (the polynomi- 
als Hj(x) multiplied by the Gaussian) peaked around a 
very high momentum ko — 37r/10 and for inertial mass 
m = 0.6. The mean value moves at the group velocity 
given by the drift coefficient v. One can notice how the 
approximation remains accurate even for small position 
spreads of few Planck lengths. For a spread a of the or- 
der of a Fermi as in a typical particle physics scenario, 
the time t needed for a significant departure would be 
comparable to many universe life-times. 



C. The relativistic regime 

In the relativistic regime for t,ra« 1 and k/m 3> 1, 
the fc-dependent Schrodinger equation (33) approaches 



the Dirac equation. The leading order and the correc- 
tions to the drift and diffusion coefficients introduced by 
the automaton evolution are 



A: 



D 



\Jk 2 + m 2 

2 

m 



1 



,2^3 



1 m 2 k 2 



6 fc 2 



-m k 



2;„2 



1 m 2 fc 4 



2fc 2 



(40) 
(41) 



The leading order in v and D correspond to the Dirac 
equation. 



D. The non relativistic regime 

In the non relativistic regime, k, m -C 1 and k/m -C 1 
the usual Schrodinger drift and diffusion coefficients are 
recovered with the following corrections 



D 



k 
m 
1 

m 



1 



1 



(42) 
(43) 



Notice that the leading terms are just the usual group- 
velocity and diffusion coefficient of the Schrodinger equa- 
tion. 



IV. TESTING THE QUANTUM AUTOMATON 

In this section we consider an elementary thought ex- 
periment for testing the granularity of the quantum au- 
tomaton based on particle flying-time, and compare it 
with the theoretical optimal in-principle testing. As we 
will see, the latter needs time duration of the order of 
an hour, whereas the former requires a time comparable 
to the age of the universe. Other experiments could be 
devised e.g. using quantum interferometry and/or ultra- 
cold atoms as in Refs. 19-22J: these will be considered 
in a forthcoming publication on the three-dimensional 
Dirac automaton. 



A. Discrimination via particle flying-time 



The fc-dependent Schrodinger equation ( 33 ) along with 



the leading terms in the relativistic and non relativistic 



regimes in subsections III C and III D provide a unique 
analytic tool for evaluating the macroscopic evolution 
of the automaton, which otherwise would not be com- 
putable in practice. As a thought experiment for an ex- 
perimental falsification of a quantum automaton evolu- 
tion at the Planck scale we consider a measurement of 
fly time of a particle. As for order of magnitude, we con- 
sider numerical values corresponding to ultra high energy 
cosmic rays (UHECR) [44]). 

Consider then a proton UHECR with m p s» 10~ 19 and 
momentum peaked around kcR. ~ 10 -8 in Planck units, 
with a spread a. We ask what is the minimal time tea, 
for observing a complete spatial separation between the 
trajectory predicted by the cellular automaton model and 
the one described by the usual Dirac equation. Thus we 
require the separation between the two trajectories to 
be greater than a = a^ 1 the initial proton's width in 
the position space. Notice that UHECR still belong to 
the relativistic regime m p , ken^ 1 (see Subsect. |HIC|, 



where the automaton well approximates the usual Dirac 
evolution. 

We describe the state evolution of the wavepacket of 
the proton using the fc-dependent Schrodinger equation 
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FIG. 5: (Colors online) Test of the approximation (381 of the Dirac automaton evolution Eq. (Ill in one space dimension. The 



approximation corresponds to the quantum-field asymptotic behaviour. Here the state ( 29 I of the Schwartz class as in QFT, is a 
superposition of Hermite functions (the polynomials Hj(x) multiplied by the Gaussian) peaked around momentum fco = 3n/10, 
specifically \ip(x,0)) = Ae ik ° x ]T\ gN Cje~ x /i& H j (x/2a)\+) ko where a = a' 1 = 20 is the position variance corresponding to 

momentum variance a, and the nonvanishing terms are Co = \/l/3, C2 = \/4/9, C7 = \/2/9. The automaton mass is m = 0.6. 
The momentum and mass parameters are in the Planckian ultrarelativistic regime. In the picture we show a comparison at four 
different times t = 0, t = 100, t = 20 an d t — 600 between the automaton probability distribution |?/>(:r,£)| 2 (in red) and the 
solution of the Schrodinger equation (33 1 \ip(x,t)\ 2 (in blue). The drift and diffusion coefficients are respectively v = 0.73 and 
D = 0.31. The mean value moves at the group velocity given by the drift coefficient v. The approximation remains accurate 
even for position spread & — 20 Planck lengths. According to Eq. (391 one has significant deviations for t w 7c- 3 , which is 
t — 600 in the present case. However, a reasonable spread a in a typical particle physics scenario is the Fermi length a ~ 10 20 , 



that would need a time t comparable to many universe life-times to introduce a significant error. The e error in Eq. (391 can 
be taken very small by considering na instead of a in Eq. (1291). For Gaussian states it is enough to consider 3a to get e w 10 -3 . 



(331 with initial Gaussian state. The Dirac evolution is 



very precisely given by the Schrodinger equation ( 33 1 just 
the leading-order terms in Eqs. (40) and (41), whereas 



the automaton is described in this regime by the full ex- 
pansion. The time required to have a separation a be- 
tween the automaton and the Dirac particle can be easily 
calculated using the power expansion in Eq. ( 40 1 



t fa a 



1 2\3 



m 2 k 2 (2m 2 + fc) 



(44) 



which, since it is m p /kcR <C 1 further simplifies as fol- 
lows 



ten « 6- 



>>)■?, 



(45) 



In order to be visible the separation a, the overall broad- 
ening abr(t) of the two packets must be much smaller 
than <t. Using Eq. (41) one has 



(Jbr(t) 




where Z?d = m 2 (fc 2 + m 2 ) 3 / 2 and we used m p /kcR <C 1. 
Using Eq. ([45]) we see that a ^> ov when 



10 22 Planck lengths = 10 2 fm, (46) 
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a physical result, considering a reasonable width for the 
proton packet of the order of a Fermi or bigger. With 
a = 10 2 fm the flying time request for complete separation 
between the two trajectories is 



i CK si6 x 10 60 Planck times « 10 iY s, 



17, 



(47) 



comparable with the universe age. In Subsec. |IVB| we 
will evaluate the minimal theoretical time required for 
discriminating perfectly between the automaton and the 
Dirac evolution using UHECR. 



B. Theoretically optimal discrimination 

The problem of discriminating optimally between two 
different dynamics is what is usually refereed to as "dis- 
crimination between two black boxes" . An experimental- 
ist is given a black box that can be either the automaton 
(box A) or Dirac equation (box D) with equal probability, 
and is asked to guess which box. The box discrimination 
is a two-outcome experiment, with outcomes A and D. 
The correct answer is given with some probability, and 
the challenge is to minimize the probability error (wrong 
guess) 



p e = ib(A|D)+p(D|A)] 



(48) 



where p(X|Y) is the probability of getting outcome X 
when the black box is Y. We want now to minimize p e 
over all the possible experiments. The optimal discrim- 
ination between the two black boxes is a special case of 
discrimination between quantum channels, namely the 
discrimination between the two unitary evolutions de- 
scribed by Ua and Ujj. Generally, the optimization be- 
tween quantum channels needs entangled states |45j , but 
for the case of two unitary evolutions with a single copy 
of each box, the optimal error probability p e is given by 

051 



p e = \-\\\U A -U^\. 



(49) 



However, such a probability with no restriction on the in- 
put states would be vanishingly small in the case of large 
dimension for the Hilbert space. We thus consider the 
physical bounds on the number of particles N < N and 
their momentum k < k. In such case the old Helstrom's 
result gives [4*o] 



Pe — 2 2 



sup \\U A pU{ - UupU^h, (50) 

where 7^ jy denotes the set of physically restricted states 

p E T^n iff Tr[pN- k ] = Tr[pP R ] = (51) 

where Pff is the projector on the N > iV-particlcs 
sector and Nj, is the operator that counts the num- 
ber of particles with momentum |fc| > k, i.e Nj, = 



J.^^dkip' (k)ip(k). In Appendix pi we evaluate a lower 
bound for the optimal probability of error, which is given 
by 



1 - cos 2 (g(k,m,N,t)) 



(52) 



where 

g(k, m, N, t) :— N arccos (cos(ai) — /?) 



a := max \ujy) 
fce{o,fc} 



(53) 



j3 := max 
fce{o,fc} 



and v (see Eq. (34)) and ud = k/^/k 2 + m 2 the automa- 



ton and the Dirac drift coefficients. 

A simplified version of the bound in the typical particle 
physics regime, namely k,m I is obtained by expand- 
ing in series the function g in Eq. (53) near m = k = 0. 



Truncating the expansion at the leading order and ne- 
glecting a small constant term we have 



f 



j{m,k,N,t) « -m 2 kNt 



(54) 



By putting p e = 0, corresponding to g(m, k, N, t) = n/2, 
we obtain the minimum time required for discriminating 
perfectly between the automaton and the Dirac evolution 



t mm {m, k,N) w 37r 



1 



m 2 kN 



(55) 



Notice that this is an in-principle optimal result, with- 
out any specification of the actual apparatus needed to 
achieve it. For a proton UHECR with k — ken ~ I0~ 8 
we have 



tr. 



minV'i'pi ^CR 



kcR, 1) ~ 3ttI0 46 Planck times 10 3 s, (56) 



which is many orders of magnitude smaller than the time 
in the flying-time experiment of Subsect. |IV A| 



V. CONCLUSIONS 

In this paper, based only on the parity and time- 
reversal invariance and without requiring Lorentz covari- 
ance, we have derived a QCA that describes the Dirac 
evolution. The automaton extends the Dirac field the- 
ory to the Planck and ultrarelativistic scales. The Dirac 
equation is recovered for small m and k. 

We have provided a technique to derive an analytical 
approximation of the evolution of the automaton in terms 
of a momentum-dependent Schrodinger equation. The 
approximation works very well for quantum states typical 
of QFT having limited bandwidth in momentum, and 
gives precise estimation of the automaton behaviour in 
all regimes, including ultrarelativistic and Planckian. 

We have then presented a thought experiment for test- 
ing the automaton using a fly-time discrimination be- 
tween the automaton state and the Dirac one, for the 
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typical UHECR scale. However, this leads to a time- 
duration of the experiment that is comparable to the uni- 
verse life-time. We have then considered the in-principle 
discrimination between the Dirac and the automaton evo- 
lution, optimizing over all measurement apparata and 
preparation available according to quantum theory, and 
providing the general discrimination formula. For the 
case of UHECR we got an experiment duration of the 
order of 10 3 s. 

The present preliminary analysis of the QCA exten- 
sion of QFT to the Planck scale has proved very effective 
in evaluating experimentable phenomenology, and this 
motivates pursuing investigation of more general QCAs, 
with possible correspondence to the theory of fundamen- 
tal interactions, and including a theory quantum gravity. 
This, in consideration of the simplicity of the automaton, 
of the principles at the basis of the QCA, and of all mo- 
tivations that we expressed in the introduction, namely 
the fact that the Lorentz covariance is not assumed a pri- 
ori but is recovered as emergent, the fact that the QCA 
is quantum ab-initio, the fact that it doesn't suffer vi- 
olations of causality, and is free from all the problems 
arising from the continuum. 

Next steps in the QCA research will be the Dirac Equa- 
tion in 3 + 1 dimensions [32], and the corresponding k- 
dependent Schrodinger equation. The 3 + 1 dimensional 
case is particularly relevant in view of the spin-statistics 
connection. At least for the Dirac automaton, also in 
3 + 1 dimensions the automaton is derived with minimal 
assumptions, as in the present case, and the usual Dirac 
equation with two spinors emerges for small m and k. As 
regards gauge-invariance the QCA has a very natural em- 
bedded quantum gauge-invariance, through the choice of 
a local unitary operator on each system of the quantum 
lattice, and this opens a route to the interacting QFTs. 
One of the problems that are still open is the possibility 
of achieving the automaton by qubits on the lattice, a 
problem that dates back to the famous work of Feynman 
[34) . In the present 1-dimensional case the simple Jordan- 
Wigner transform achieves the transformation from lo- 
cal qubits to nonlocal Fermi fields [3§]- The extension 
to higher dimensions of the Jordan- Wigner is straight- 
forward, however, it is still not clear if it is possible to 
map local field interactions to local interactions between 
qubits. It seems that a possible solution requires auxil- 
iary fields of the Majorana kind [251 S3 EE] ■ The pos- 
sible solution of the Feynman problem may shade light 
on the intimate difference between the Fermionic and the 
Bosonic natures. 



Appendix A: Derivation of Eq. (39) 



Here we evaluate the overlap between the ex- 
act automaton evolution \^p(t)} and the fc-dependent 
Schrodinger approximation \ip(t)) 



\(mm))\ = 



> 



2lT 



2vr 



* *e+ l << )fe3 +°( fc4 »<| ff (fc,0)| 2 

27T 



dke~ l( < k3+ ° ik4))t \g(kM 2 



> 



ka — a 



djfce- i( < )fc8+o(fc4))t | S (Jfe,0)| 2 

\k— k()\>a 



> 



/ dk\g(k,0)\ 2 -0(a 5 )t 

3+ <nf~5\ 



- e > 



> 1- e- 1( j A t-0{cj b )t 



with the constant 7 = J k °_ a dk \g(k, 



Appendix B: Proof of the bound (52 1 



In this appendix we detail the proof of the bound ( 52 ) 
in Section flV B| which provides the probability of optimal 
error probability in discriminating the Dirac automaton 
and the usual Dirac evolution. The discrimination ex- 
periment can have a generic duration t and the unitary 
operators to be discriminated are explicitly given by 



U'(fc) = cxp(-iH(k)t) 



■1 s 

-10- 



COS(LOt) + % ^— L 0' 

—to — * — 1 cos(w£ — 1 — - — ~( 



(Bl) 



U^(fc) = exp(-iH D (fc)t) 



cos(Ai) + i 



sin(At) 



cos(Ai) 



q(At) > 
■ sin(At) j t 



(B2) 



as can be easily verified by direct computation using the 
Hamiltonians in Eqs. (14) and (20). The proof of the 



bound ( 52 ) goes throug 



1 the following three Lemmas. 



Lemma 1 LetU^^k) andU*^(fc) be defined according to 
Eqs. prllpLl and let us define V(k,t) = U^U^O) 



Let e l ~ m fl (, c an eigenvalue ofV(k,t). Then the fol- 
lowing bound holds: 
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cos(^(fc, to, t)) > cos(at) — f3 



where 



(B3) 
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a(k, to) 



P{k,m) := - [1-vvn- <J (1 - v 2 )(l - v& 



(B4) 



11 



Proof. Since both U^(fc) and U tf (fc) are SU(2) ma- 
trices, we have that V(fc,i) is an SU(2) matrix and its 
eigenvalues must be of the form e v( fc . m >*) and e -v( fe > m . 4 ). 
This implies the equality cos(/i(fc, to, t)) = |Tr[V(fc,t)] 
which by direct computation gives 



P 



P 



cos(fj,(k,m,t)) =1 I cos(ai) H — cos(7i) (B5) 



which is satisfied if 

w(k, to) := to 2 cos(fc) + sin 2 (fc) - nu^ > . (B16) 
It is easy to see that, for any to g [0,1], we have 
(d^wik, to)) = for i = 0, 1, while d[ 2) /(fc, rn) > 

V / k— 

for any fc € [0,7r/2], which gives the monotonicity of /?. 



where a and /3 are defined accordingly with Eq. (B4) 
and 7 := lu + cjrj- Finally, from Eq. (B5| one has the 
bound cos(/i(fc, to, t)) > cos(at) — /? ■ 

The second Lemma shows the monotonicity of the two 
functions a,/3 in Lemma [l] 

Lemma 2 Let a(k,m) and f3(k,m) be defined as in Eq. 
(B4| and < fc < tt Then we have 



Lemma 3 Let < fc < tt, N be a positive 
number, and a, ft be defined as in Eq. (B6). 
1 — cos(^y ) and t < /(fc, m, N) where 



f(k,m,N) := 



arccos(cos (^) + P) 



then 



a := max_ |a| 

ke[-k,k] 

(3 := max \f3\ 

ke[-k,k] 



max \a\ 

ke{0,k} 

max |/3| 
fce{o,fe} 



Nfi(k, to, t)) < g(k, to, iV, t) < 



Vto e [0, 1] . (B6) w /j ere g fa TOj TV, t) := N arccos (cos(at) - /?) . 



Proof. Since both u> and wd are even function of fc, from 
Eq. (B4) we have that also a and /3 are even function 



of k. For this reason we can restrict to k £ [0,fc]. The 
equality ( B6 1 can be proved by showing that a and /3 are 



nondecreasing function of k for fc € [0, A;]. 

Since dkCt = vr> — v, clearly v 2 ^ — v 2 > for k £ [0, tt) 
implies dkot > in the same interval. By direct compu- 
tation one can verify that 



x(k, to) 



y(k,m) 



y(k, to) 

x(k, to) := k 2 - sin 2 (fc)(l - to 2 ) 
(fc 2 + TO 2 )(sin 2 (fc) + m 2 cos 2 (fc)) . 



(B7) 

(B8) 
(B9) 



Clearly we have y(fc, to) > and since k > sin(fc) for 
< k < tt, the thesis is proved. 

Again the monotonicity of f3 for k £ [0, w) follows from 
dkP > in the same interval. By elementary computa- 
tion we have 

dkP = x(k, m)y{k, m)z(k, to) (B10) 

2 

TO 

x(k,m):= ^— (Bll) 

wd sin (u) 

y(k,m) := (nsm(fc) - fc) (B12) 

ncos(fc) 1 

z(k,m) := ^ B13 

sin (w) 

Clearly x(k,m)y(k,m) < for fc G [0,tt) and we just 
have to verify that z(k,m) < in that interval, namely 



to 2 cos 2 (fc) + sin 2 (fc) — n cos(fc) > . 



(B14) 



The last equation is clearly satisfied for fc £ [tt/2,it] 
therefore we restrict to fc £ [0, tt/2}. This allows to divide 
the left side of Eq. (B14) by cos(fc) achieving 



to 2 cos(fc) 



sin 2 (fc) 
cos(fc) 



— nuj D > 



(B15) 



Proof. The conditions t < f(k,m,N) and j3 

< at < arccos(cos (^y) ^ = ^ 
=>■ 1 > cos(ai) - /3 > cos ( =>■ 



integer 
IfP< 



(B17) 
(B18) 

< 1 - 



cos(ai) — P > cos(at) — P > cos (^y 7 ) 



(B19) 

By exploiting the bound (|B3|) into Eq. (|B19|) we have 



cos(/i(fc, to, t)) > cos(ai) — P > cos 



Nfj,(k, m,t) < N arccos (cos(at) — P) 



(B20) 



< 



We are now ready to prove the bound ( 52 ) 



Proposition 1 Let U f and be the unitary evolution 
given by the Dirac QCA and by the Dirac equation re- 
spectively. If the hypothesis of Lemma [3| hold we have 



SUp WiU'pU^ -U^pU^Wt < ^l-CO S 2(« ? (fc,TO,iV,<)). 

(B21) 

Proof. First we notice that thanks to the convex- 
ity of the trace distance we can without loss of gen- 
erality consider p to be pure. If p is a pure state 



\x)(x\ the trace distance becomes yl— KxI^'^d \x)\ 2 = 
y/l — \(x\V(t)\x)\ 2 - If we expand \x) on a basis of eigen- 
states of V, i.e. \x) = Z)iv,k,s \/PN,k, s \N,k,s), we have 



\(x\v(t)\ x }\ 



N 

PN,k, s exp I i^2sjn(kj,m,t) 

.V.k.s \ j=0 



E 



> 



> 



E 



N 



PJV,k,s COS ■ ^ o 3 
W,k,s \j=0 



Sj(j.(kj,m, t) 



(B22) 
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By exploiting the bound (B18) into Eq. (B22) we have 

2 



N 



^ P7v,k, s cos y^Sj/x(fcj,m,f) 
which finally implies 



> cos 2 (g(k, m, N, t)) 



Vl-\(x\V(t)\ X )\ 2 < Jl - cos 2 ( 3 (fc, m, #, *)) . 



Inserting the bound ( B21 ) into Eq. ( 50 ) we finally have 
the bound ( 52 ) . 
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